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Abstract 

The Ginzburg-Landau theory of a trapped Fermi gas with a BEC-BCS crossover is derived by the 
path-integral method. In addition to the standard Ginzburg-Landau equation, a second equation 
describing the total atom density is obtained. These two coupled equations are necessary to de- 
scribe both homogeneous and inhomogeneous systems. The Ginzburg-Landau theory is valid near 
the transition temperature Tc on both sides of the crossover. In the weakly-interacting BEC region, 
it is also accurate at zero temperature where the Ginzburg-Landau equation can be mapped onto 
the Gross-Pitaevskii (GP) equation. The applicability of CP equation at finite temperature is dis- 
cussed. On the BEC side, the fluctuation of the order parameter is studied and the renormalization 
to the molecule coupling constant is obtained. 
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I. INTRODUCTION 



In a Fermi gas with a BEC-BCS crossover, the scattering length can be tuned by the 
technique of Feshbach resonance |l| . The effective interaction between atoms is proportional 
to the scattering length. A dilute Fermi gas with a negative scattering length is in a BCS 
pairing state below a critical temperature 2, s], ^, Q], very similar to BCS superconductors. 
A dilute Fermi gas with a positive scattering length is in a BEC state of diatomic molecules 
below another critical temperature Although the scattering length is divergent at 

the resonance, the system evolves smoothly between the BCS state and the molecular BEC 
state across the resonance at low temperatures. The observation of the BEC-BCS crossover 
8| provided a new platform to study strong-correlation effects in fermionic systems. 

The BEC-BCS crossover can be qualitatively understood in the BCS-type mean-field 
theory jo, [lo|. In this theory, as the interaction changes sign from attractive to repulsive 
across the resonance, the pair size of Cooper pairs decreases, and eventually these atom 
pairs become diatomic molecules. Although the mean-field theory offers the correct physi- 
cal picture, it overestimates the critical temperature and the molecule-molecule scattering 
length in the weakly-interacting BEC limit. Nozieres and Schmitt-Rink (NSR) ll| found 
that fluctuation effects have to be considered to get the correct critical temperature. The 
total fermion density includes not only the mean-field fermion density but also the density of 
thermal molecules, which provides the important relation between the density and the chem- 
ical potential. The NSR theory is essentially equivalent to treating Gaussian fluctuations 
■n the Ginzburg-Landau theory Q. The NSR theory was also apphed at zero temperature 
and the molecule-molecule scattering length was found in good agreement |13l . |14I | with the 



few-body calculation []J 

The purpose of this paper is to construct the Ginzburg-Landau theory to describe 
the BEC-BCS crossover in a trapped Fermi gas. Compared to microscopic theories, the 
Ginzburg-Landau theory has potential advantages of requiring less computation and be- 
ing easier to be applied to inhomogeneous cases such as trapped systems. In the weakly- 
interacting BEC region, the Ginzburg-Landau equation was shown to be equivalent to the 
Gross-Pitaevskii (CP) equation at zero temperature [l^. In the unitary region, a mod- 



17|, vortex [l8|, 



ified Ginzburg-Landau theory was developed to describe the phase slip 
and vortex lattices 1^. However there still lacks a complete Ginzburg-Landau description 
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of the whole BEC-BCS crossover. In the following, we first derive the Ginzburg-Landau 
theory of a trapped Fermi gas by the functional-integral method, and obtain the Ginzburg- 
Landau equation and the equation for the fermion density. The density equation is impor- 
tant for providing the density profile of the BEC-BCS crossover in both the inhomogeneous 
and homogeneous cases. Then we concentrate on weakly-interacting BEG limit, study the 
Ginzburg-Landau equation at both zero and finite temperatures, and consider effects due to 
fluctuations of the order parameter. The conclusion is given in the end. 



II. GINZBURG-LANDAU THEORY OF A TRAPPED FERMI GAS 



A Fermi gas with a wide Feshbach resonance can be effectively described by a single- 
channel model, while for the narrow resonance case a two-channel model is more accurate 
Q. In this pape. we cos.de. only the wide —ce case which the s,ngle-cha,.el 
Hamiltonian density is given by 



X 



2m 



+ Viv) 



(l)a{x) + g(f)\{x)(f)\{x)(f)i{x)(f)^{x), 



(1) 



where x = (r, r) is the coordinate in space and time, (f)a{x) is the fleld operator of atoms with 
spin-component a, m is the mass of a Fermi atom, ^(r) is the external trapping potential, 
and the coupling constant is given hj g = Anfi^as/m. In the following we consider only the 
spin-balanced case where the densities of spin-up and spin-down atoms are the same. 
The grand partition function Z can be written in functional-integral formalism as 



(2) 



with the action given by 



(3) 



where < r < 1/(A;bT), and fi is the chemical potential. 

The interaction term in Eq. ([T]) can be decoupled by introducing an auxiliary fleld 
A{x) and applying the Hubbard-Stratanovich transformation 12|]. After integrating out the 
fermion fleld (i>(x), we obtain 



^A*^Aexp[5A] 



(4) 



where the action S'a is given in terms of the auxihary field A, 

5A = lndetM + - [ d^x\A(x)\\ 
9 J 



M 



-dr - H{y) -A(x) 
-A*(x) -dr + Hir) 



and 



2m 



Hir) 



V{r) - /i. 



The action S'a can be further separated into two parts, 

•S'a = Seff + So, 



(5) 
(6) 

(7) 
(8) 



where 5*0 = In det Mq is independent of A with 

-dr - H(r) 

Mn 







-dr + H{r 

and the second part Sejf vanishes when A = 0, 



(9) 



with u = Mo - M 



d\|A(x)|^ + Trln[I-Gu], 
I as the identity matrix, and G = Mq ^ 



(10) 



G+ 
L G_ 



A (a;) 
LA*(x) 

being the Green's function of a noninteracting Fermi gas. 

Next we expand the second term in Seff to the fourth order in A. This approximation 
holds only when A is small compared to the Fermi energy, which is true when \T—Tc\/Tc <^ 1 
or the system is in the weakly-interacting BEG regime. After the expansion, we obtain 

..„.i/d%|AWP + /d^.d*.,Q(...,)A-WA(.,) 

1 r ^ 

~9 / n^'^^^^^i'-'^4)A*(a:i)A(a:2)A*(a:3)A(x4), 



(11) 



where 



Q{xi,X2) = -G+{xi,X2)G-{x2,Xi), 

...,X4) = G+{xi,X2)G-{x2,Xs)G+{x3,X4,)G-{x4,Xi). 



(12) 



Since we are interested in low-energy and long-wavelength properties of the system, we 
apply gradient expansion in Eq. f|TT]) and obtain 



S, 



dA*{x)drA{x) + cA*{x)^-^A(x) 



Am 



where 



+ {a + \)\A{x)\'-\h\A{x)\' 
9 ^ 



d^x'Q{x-x'/2,x + x'/2), 
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Y\_d'^XiR{x,Xi,X2,X3), 
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(13) 



^ I d'x'^—Q{x-x'/2,x + x'/2), 



(14) 



and x' = (r', t'). The coefficient d of the time-derivative term is defined as the coefficient of 



the linear term in the expansion of Q in the frequency space in the zero-frequency limit 

' - 1 



d = lim / d^x' 



'-Q{x-x'/2,x + x'/2). 



(15) 



The equation of motion of the order parameter can be obtained by taking 6Seff /SA*{x) = 0, 
which yields 

(16) 



(ddr + c + a + -)A(x) - b\A(x)\^A(x) = 0. 
4m g 



For an arbitrary trap potential V{r), it is difficult to obtain the exact analytical expression 
of the Green's function G. Here we consider only the case where the trap length is much 
larger than the inter-particle distance, and the local-density approximation (LDA) can be 
applied. 



G{xi,X2)^G^"\x,,X2). 



(17) 



Here G^o) 







is the Green's function of a homogeneous noninteracting Fermi 





gas with the chemical potential fi' = fi — V[{ri + r2)/2] 



-dr-T + fi' 

-dr + f-fi' 



G^^\x-x') = Six-x% 



(18) 
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where T = —h?'V'^/{2m). In LDA, the coefficients in Eq. fll4l) are given by 



d^A; tanh[ek/(2A;BT)] 


(27r)3 


2ek 




'tanh[ek/(2A;BT)] 


(27r)3 


[ K 




"tanh[ek/(2A;Br)] 


(27r)3 




d^k tanh[ek/(2A:BT)] 


(27r)3 
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SkBTel 



S/csTek 

^ = / — , (19) 



where ek = h'^k'^/{2m) — fi'. It is important to note that in obtaining the coefficients given 
by Eq. (fT9l) LDA is applied to the Green's function of a noninteracting Fermi gas G^'^'', 
which is different from applying LDA directly to the broken symmetry state. The latter 
case requires that the gap has to be much larger than the trap frequency, which does not 
hold near the critical temperature or in the weakly-interacting BEG regime. 

In the unitary and BGS regime where there is a clear Fermi surface, the Fermion energy 
ek has zero points, and there are divergences in integrands on r.h.s. of Eq. (fT9|) at zero 
temperature, suggesting that this formalism does not work at zero temperature. At finite 
temperatures, the integrand on r.h.s. of the equation for the coefficient d also diverges. 
A more careful treatment of this coefficient leads to damping in the dynamics of the or- 



der parameter 



211 ]. In the following, we will not study the dynamics or zero-temperature 



properties in this regime so these issues do not occur. 

At the stationary state, the function A(r) satisfy the saddle-point condition 



0, (20) 



5A{r) 

which leads to the Ginzburg-Landau equation 

- c|^V'A(r) - (a + ^)A(r) + 6|A(r)pA(r) = 0. (21) 

The Ginzburg-Landau equation can be rewritten in terms of a wave-function tp{r) defined 
by tpir) = ^/cA{r), 

- + a^(r) + /5|^(r) I V(r) = 0, (22) 



where 



a = — [a -\ — 
c 9' 



P = (23) 



C2 



and the gradient terms of c are ignored due to its small spatial variation within LDA. In 
Eq. fl23l) . the coupling constant is given by 



^ ' Airh^a, j (27r)3 2e0' ^^^^ 

where = /i^/i;^/(2m), and the second r.h.s. term is a counter term in the particle-particle 
channel. 

In the simple homogeneous case, the critical temperature Tc is determined from the 
equation 

a = 0, (25) 

which is exactly the Tc equation in the mean-field BCS theory, 

d^A; r tanh[ek/(2A;BTe)] 1 ) m 



d^A; 1 



(27r)3\ 2ek 26°/ AixK^as' ^^^^ 

Above Tc, the Ginzburg-Landau equation does not have a nontrivial solution. Close to T^,, 
to the first order of T — T^, the coefficient a is approximately given by 

a« A(TJ(T-Te), (27) 

where 

d^k 1 , 2r 1 



m) = 7TT-TI7-r^sech^ 



(27r)3 4cA;BT2 '2kBTc'' 
In the trapped case, the critical temperature is approximately determined by 

a(ro) = 0, (28) 

where ro is the place with the highest fermion density. 

The Ginzburg-Landau equation determines the distribution of the order parameter in the 
stationary state. To obtain the density distribution, we need to study the thermodynamic 
potential Q, given by 

n = -ksTlnZ = Qf + ns, (29) 
where Qf is thermodynamic potential of a noninteracting Fermi gas, in LDA it is given by 

/f d^k 
dV J [ek + ksT In rik] , (30) 

and nk = 1/{1 -|- exp[ek/(/i:B^)]} is the Fermi distribution function. The term Qg is the 
contribution to thermodynamic potential by the order parameter. If we consider only the 



stationary state and ignore fluctuations of the order parameter in the mean-fleld approxi- 
mation, Qs is approximately given by 



eff 



(31) 



From the relation = —dQ/dfi, the atom density in the mean-fleld approximation can be 
obtained, 

n{r) ^nf{r)+nf\r), (32) 

where Uf associated with Qf is the density of a homogeneous noninteracting Fermi gas with 
the chemical potential /i' = /i — V^(r) , 

/d^k 
(^'^k(r), (33) 

and nf^ associated with l^i^-* is the density due to the order parameter, 

nf(r)^|^|A(r)r = 2|^(r)|l (34) 

Note that in Eq. (134|) the quartic term in ip{r) is ignored because it is much smaller than the 
quadratic term. The density equation fl32|) indicates that the total density can be separated 
into two parts, i.e. Uf from the normal state and n^^^ from the superfluid order parameter, 
which is consistent with the two-fluid model of a superfluid at flnite temperatures. The 
superfluid atom pairs are described by the wave-function 4j{r). In a trapped system, usually 
the total number of atoms is given, 

N = I dVn(r), (35) 



from which the chemical potential fi can be solved. The Ginzburg-Landau equation (!22|l 
and the density equation fl32|) provide a complete mean-fleld phenomenological description 
of a trapped superfluid Fermi gas. 



III. GINZBURG-LANDAU EQUATION IN THE WEAKLY INTERACTING BEG 
REGIME 

In the weakly interacting BEG regime, the Fermi gas is dilute, na^ -C 1, and the order 
parameter A is much less than the binding energy of a diatomic molecule given by eo = 
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h'^/{'ma1). In this regime, since the chemical potential /i is negative, there is no Fermi 
surface. The coefficients in the Ginzburg-Landau equation given by Eq. f|T9l) are well 
defined even at zero temperature. Therefore the Ginzburg-Landau theory can be applied 
from zero temperature to near T^. In this section we study the Ginzburg-Landau equation 
in this regime. 

At zero temperature, from Eq. (fTOjl and (l23l) . the coefficients in the Ginzburg-Landau 
equation are given by 



ao = A^/\f7\{^/\fJ'\ - A/eo/2), 

^° ~ 3/2 /r^' 

do = Co. (36) 

Since the density tiq = |«o//5o| is much smaller than l/a^, we obtain /i' = — (eo/2)[l + 
O(raoag)], and approximately 

ao ^ -(2/i' + eo), 

A « (37) 

2m 

With the time-dependent term, the Ginzburg-Landau equation can be written as 

dM^) - ^H^) + [2^(r) - f^bW^) + = 0, (38) 

where fib = '^fJ' + ^o is the chemical potential of molecules. This time-dependent Ginzburg- 
Landau equation is identical to the Gross-Pitaevskii equation of molecular BEG if the imag- 
inary time T is analytically continued to the real time t, t = iht, with = 2m identified 
as the molecule mass and ab = 2as identified as the scattering length between molecules. 
The density equation ( l32l) in this regime is trivial, n = 2|'?/;p, which means all the atoms 
are paired into condensed molecules at zero temperature. However, the molecule scattering 
length ttb extracted from Eq. (!38|) is = 2as, contradicting to the result ~ O.Ga^ from the 
few-body calculation [l^. This discrepancy is due to the fact that we have not considered 
the fiuctuation effect which are discussed in the next section. 

In the weakly interacting regime, the BEG transition temperature Tc is much smaller than 
the molecule binding energy eo, ksTc ^ eo. As a result, at any finite temperature below 
or near T^, the coefficients in the time-dependent Ginzburg-Landau equation are almost 
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the same as those at zero temperature given by Eq. (138!) . except the molecular chemical 
potential fib now varying with temperature. Therefore the time-dependent Ginzburg-Landau 
equation (1381) remains valid at finite temperatures below or near T^.. 



IV. FLUCTUATION EFFECTS 

In this section, we consider the effect due to the fluctuation of the order parameter which 
is ignored in our derivation so far. The effective action in terms of the pair wave-function 
■ip{x) = a/cA(x) is given by 



Seff = I d X 



Kllj*(x)dr^lj(x) + tlj*(x)-^V'^tlj(x) 

Am 



(39) 



where k = d/c. For simplicity we consider only the homogeneous case. In the following, we 
concentrate on the weakly- interacting BEG regime near or below Tc, where k k, 1^ a ^ ~f^b, 
and (3 ~ 87r^^as/(2m). In this regime, the effective action given by Eq. (139|) is the same 
action of a Bose gas with the boson mass given by 2m and a scattering length given by 2as. 

The effective action given by Eq. (l39l) provides a contribution to the thermodynamic 
potential given by 

ns = -ksTln j ^i)*S)i) exp[5e//], (40) 
from which we can obtain its contribution to the density n .o by taking —dQs/dfi. For a 



weakly-interacting Bose gas, Bogoliubov's theory is accurate 



in which Ug is given by 



o,2 , f r(efek-tt) Ebk . . , . 



where e^k = fv^k"^ / (4m) is the kinetic energy of a molecule, and -E^k = a/ ^hki^bv. — 2a) is the 
excitation energy of the molecular quasi-particle. The order parameter ifjQ is equivalent to 



the expectation value of the Bose fleld-operator, = ^J—a/ 13, where its phase is chosen so 
is positive for simplicity. 

In Bogoliubov's theory of a dilute Bose gas, there is an ultra-violet divergence appearing 
in the calculation of the thermodynamical potential VLg and the ground state energy, which 
is removed by the renormalization of the coupling constant 

_i 2m f d?k 1 

9b - 



AnK^ab J {2n)^ 2ebi,' 
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where the second r.h.s. term is a counter term only appearing in the calculation involving 
the particle-particle channel. However in the effective action Eq. fl39|) . the constant (3 is 
not renormalized so far, which would result in a divergent term in VL^ and the ground state 



energy 



231] given by 



k 



with the contribution to density given by 



[ d?k I 

where dajd^i ~ —2. This ultra-violet divergence is an unphysical result, due to the invalidity 
of the gradient expansion used in derivation the Ginzburg-Landau action given by Eq. ( |T3|) 
at short distances. Fluctuations inside molecules become important at short distances, 
which is beyond the description of the gradient expansion. The energy scale at which 
the gradient expansion is invalid is approximately given by the molecule binding energy 
^0 = ^ I ("^<2^), which is equal to the molecular kinetic energy at wavevector k = 2/ as- Thus 
a straightforward renormalization method is to put a cutoff A = 2/as in the /c— integrals of 
Eq. (02]) and (03]), which yields 

After considering the renormalization due to fluctuations, we obtain the density equation. 



" - + + 2(1. + / |||(-^„th(^)-tanh(5g^)|^(45) 

At zero temperature, the density equation is simply given by 

n = 2(l + -)V'o+ / — ^ 1] 

= 2(1 + ^)^o^ + (46) 

The first r.h.s term in Eq. ( l46i) implies that the condensate density is rio = 2(l + 8/7r)V'o not 
simply 2iIjq. The second r.h.s term is proportional to tpQ, which comes from the quantum 
depletion of molecules. Thus in term of the true molecular condensate wave-function 

^P{x) = v/(l + 8/7r)^(x), (47) 

the Ginzburg-Landau equation given by Eq. fl38|l should be rewritten as 
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FIG. 1: The molecule scattering length ai, vs the cutoff A. The dotted lines are A = 2 /a., and 



fflfe = 0.56os which is very close to the result at 0.6as from the few-body calculation 
A = 0, the mean-field result a{, = 2as is recovered. 



At 



dMx) - -^^(x) + [2V(r) - fi,Mx) + ^— = 0, (48) 
where a;, is the scattering length of molecules after considering the renormalization, 

ab = ~ 0.56a„ (49) 

1 + »/7r 



very close to the result ab ~ CGa^ from the few-body calculation pj 

Although the renormalization to the mean-field molecule scattering length given by Eq. 
(jlHj) was obtained below Tc, it is valid above Tc as well, which can be understood in the 
vacuum renormalization of the molecule coupling constant in the T-matrix approximation, 

9b^ = 9m y ^^G'f,(k,^)Gf,(-k,-^) 

where Qbo = Sir hj^ag/ {2771) is the molecule coupling constant in the mean-field approximation, 
gb = 47r/i^afe/(2m) is the renormalized coupling constant, and 

G,{k,u;) = ^— - 

— Cfek + M 
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is the Green's function of a molecule in vacuum. A cutoff A = 2 /as should be put in the 
/c— integral of Eq. fl50l) due to the same reason as stated above that this effective description 
cease to be accurate when e^k > Cq, leading to the same renormahzation, 

= l+4a.A/7r ' ^^^^ 

same as Eq. (H9l) at A = 2 /as- The molecule scattering length as a function of the 
cutoff A is shown in Fig. [H It should be emphasized that a more precise cutoff than 2/a<j is 
necessary for any further numerical comparison with the few-body result, which requires un- 
derstanding of high-energy processes beyond the Ginzburg-Landau description. The vacuum 
renormahzation of the molecule coupling constant can also be applied to trapped systems 
within LDA as long as the molecule binding energy is much bigger than trap frequencies. 
At Tc, a = 0, the density is given by 

_ /• d^A; 1 1 

J (2^^exp[66k/(A;BTe)] - 1 + exp[6k/(fcBr,)] + 1^' ^ ^ 

indicating that there are only thermally-excited atoms and molecules. Since /c^Tc <^ eo, the 
number of thermally-excited atoms are negligiblly small, and almost all the particles are 
thermal molecules, 

"""^ J (27r)3exp[6,k/(A;BTe)]-l- ^ ^ 

Generally at finite temperature below T^, the density can be separated into the superfiuid 
density Usj and the normal density n„, n = Ugf + nn- The superfiuid density is given by the 
condensate density and the quantum depletion. 



n., = 2^o^ 1 + , (54) 




consistent with traditional theories of a dilute Bose gas [22| . The normal density is given by 



the total density of thermal atoms and molecules. Since the thermal atoms are negligible, 
the normal density is approximately given by the density of thermal molecules, 

d^k (cbk -a) 1 



(55) 



(27r)3 Efek exp[E,^/{kBT)]-l 

Compared with the mean-field result in Eq 0321) , both the superfiuid density and the normal 
density are renormalized. In the weakly-interacting BEG limit, as shown in Fig. [21 the 
renormahzation to the density is quite strong. At zero temperature, the mean-field density 
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TIT 
c 

FIG. 2: Density ratios of a homogeneous Fermi gas in the weakly-interacting BEC Umit as functions 
of temperature. The sohd hne is the mean- field density given by Eq. (|32p divided by the total 
density. For comparison, the dashed line is the ratio of the superfluid density to the total density. 

account for only about 28% of the total density, while the rest density is due to fluctuation 
contribution. At Tc, almost all the density is due to fluctuation contribution. 

Away from the weakly-interacting BEC limit, the fluctuation effect is more difficult to 
deal with in the Ginzburg-Landau theory. The coefficients in the Ginzburg-Landau equation 
are more complicated than those given by Eq. (136!) . and fluctuations of the order parameter 
cannot be simply treated by the Bogoliubov's theory. Moreover the wavevector cutoff 2/ as 
vanishes in the unitary region where a more subtle renormalization scheme is required. These 
problems will be explored in our future work. The situation is simpler again on the other 
side, in the weakly-interacting BCS limit, where the fluctuation of the order parameter is 
strongly damped. In this limit, the mean-field theory is accurate and the fluctuation of the 
order parameter is less important. 

V. CONCLUSION 

In conclusion, we have derived the Ginzburg-Landau theory of a trapped Fermi gas with 
a BEC-BCS crossover. Two equations including the standard Ginzburg-Landau equation 
and the density equation are obtained to describe the order parameter distribution and 
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the density profile. In the weakly-interacting BEC limit, the Ginzburg-Landau equation is 
equivalent to the Gross-Pitaevskii (CP) equation. The fluctuation of the order parameter is 
strong in this limit, which can be treated by the Bogoliubov's theory. Compared with mean- 
field results, both the density and the molecule-molecule scattering length are renormalized, 
in agreement with the few-body and NSR theories. This work is supported by NSFC under 
Grant No. 10674007, and by Chinese MOST under grant number 2006CB921402. 
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